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Abstract 

The elasticity and piezoelectricity of zinc oxide (ZnO) crystals and single layers are investigated 
from the first-principles calculations. It is found that a ZnO thin film less than three Zn-0 layers 
prefers a planar graphite-like structure to the wurtzite structure. ZnO single layers are much more 
flexible than graphite single layers in the elasticity and stronger than boron nitride single layers 
in the piezoelectricity. Single-walled ZnO nanotubes (SWZONTs) can exist in principle because 
of their negative binding energy. The piezoelectricity of SWZONTs depends on their chirality. 
For most ZnO nanotubes except the zigzag type, twists around the tube axis will induce axial 
polarizations. A possible scheme is proposed to achieve the SWZONTs from the solid-vapor phase 
process with carbon nanotubes as templates. 

PACS numbers: 62.25.+g, 62.20.Dc, 77.65.-j 
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I. INTRODUCTION 



Zinc oxide (ZnO) materials have attracted extensive attention for half a century be- 
cause of their excellent performance in optics, electronics and photoelectronics.— They are 
important for low cost productions of green, blue- ultraviolet , and white light-emitting de- 
vices due to their wide band gap (~3.37 eV) and large exciton binding energy (~60 meV). 
They can also be used as sensors and transducers owing to their strong piezoelectricity. 
Recently, many ZnO nanostructures have been synthesized in experiments.— Among them, 
the quasi-one-dimensional structures have become the leading research objects because of 
their novel chemical, electrical, mechanical and optical properties.— The ultraviolet lasing in 
ZnO nanowires^ and strong photoluminescence in ZnO nanorods^ have been demonstrated. 
Nanobelts, nanorings and nanohelixes are also synthesized Ji?!^ 1 ? which may be useful for 
field-effect transistors.— 

Since the discovery of carbon nanotubes in 1991,— the synthesis of tubular nanostruc- 
tures has raised worldwide interest. A lot of inorganic nanotubes, such as ZrS2, NbSe2, 
Si0 2 , Ti0 2 , BN nanotubes, etc. are achieved by several groups.— Researchers have also 
tried to fabricate ZnO nanotubes through various methods including thermal reduction,— 
vapor phase growth , 15 i 16 i 17 hydrothermal growth,— vapor-solid process,— sol-gel template 
process,— plasma-assisted molecular beam epitaxy,—^ pyrolysis of zinc acetylacetonate,— 
and Zn(NH 3 )4 + precursor thermal decomposition.— All ZnO nanotubes thus obtained have 
the wurtzite structure with diameter of 30-450 nm and thickness of 4-100 nm. Much effort 
has been made to realize much thinner and smaller ZnO nanotubes. A natural question 
is: Can we manufacture single- walled ZnO nanotubes (SWZONTs)? In this paper, We will 
answer this question based on the first-principles calculations (the ABINIT packag and 
the experimental method a 15 i 16 i 17 i 19 for synthesizing ZnO nanobelts and nanotubes with the 
wurtzite structure. Additionally, we will investigate the elasticity and piezoelectricity of ZnO 
crystals and single layers within the framework of density-functional theory (DFT) 2S and 
density-functional perturbation theory (DFPT) . 27 ' 28 ' 29 The piezoelectricity of single- walled 
ZnO nanotubes is also derived from the piezoelectricity of ZnO single layers. 

The rest of this paper is organized as follows: In Sec. m we present the results of elastic 
constants and piezoelectric coefficients of ZnO crystals. We also compare them with the 
experimental results and previous theoretical results in literatures. The main aim of this 
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section is to test whether our computational procedure can validly describe the physical 
properties of the ZnO system or not. In Sec. IIII1 we optimize a Zn-0 single layer cut from 
the wurtzite crystal and arrive at a planar graphite-like structure. Comparing the binding 
energy between the crystal and the planar single layer, we deduce that a ZnO thin film less 
than three Zn-0 layers prefers a planar graphite-like structure to the wurtzite structure. The 
elastic and piezoelectric constants of ZnO single layers are also calculated. In Sec. IIVI we 
calculate and optimize different SWZONTs. The SWZONTs are shown to exist in principle 
because of their negative binding energy. Their piezoelectricity depends on the chirality, 
and for most of them except the zigzag type, twists around the tube axis will induce axial 
polarizations. We also propose a possible scheme to realize the SWZONTs from the solid- 
vapor phase process with carbon nanotubes as templates. In Sec. |VJ we give a summary 
and propose some potential applications of SWZONTs. 

II. ZINC OXIDE CRYSTALS 

The ground state of a ZnO crystal has a wurtzite structure (space group P6 3 mc) with 
two zinc and two oxygen atoms per unit cell. The lattice constants and internal parameter 
obtained from experiments £1X6 clS follows: a = b = 3.250A, c = 5.204A, and u = 0.382, 
respectively.^ 

We optimize the structure parameters and calculate the elastic and piezoelectric constants 
of ZnO crystals. The calculations are carried by taking Troullier-Martins pseudopotentials,— 
plane-wave energy cutoff (ecut) of 60 Hartree, and 6x6x3 Monkhorst-Pack k-points^ in 
Brillouin-zone. The exchange-correlation energy are treated within the local-density approx- 
imation in the Ceperley- Alder form^ with the Perdew-Wang parametrization.— 

Using the experimental parameters to construct an initial unit cell, we have optimized 
the structure parameters and summarized them in Table HI Results by other groups are also 
listed for comparison. It is easy to see that our DFT results are quite close to those of Wu 
et al. and agree well with the experimental results (errors in 2%). 

Adopting the optimized structure, we can calculate the elastic constants and piezoelec- 
tric coefficients of ZnO crystals. The elastic constants reflect the stress-strain relation of 
materials. In terms of the symmetry of ZnO crystals (wurtzite structure), this relation can 
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be expressed in the matrix form: 
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where <Tj and q (i=l,- • -,6) represent the stresses and strains, respectively. There are only 
five independent elastic constants: Cn, Ci 2 , Ci 3 , c 33 , and c 44 . Their values in the present 
work are listed in Table ITT1 

Similarly, the piezoelectricity can be expressed in the matrix form:— 
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where Pi, P 2 , and P 3 are three polarization components along a direction, the direction 
perpendicular to a and c, and c direction, respectively. P 3 is the spontaneous polarization 
along c direction. There are only 3 independent piezoelectric coefficients: e 3 i, e 33 , and e^. 
It is necessary to point out that here the definition of piezoelectricity is different from that 
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in Ref. |37| where the piezoelectricity reflects the relation between polarizations and stresses. 
The piezoelectric coefficients obtained by the present work are listed in Table 11111 

Good agreements have been achieved between the present results and those in literatures 
(see Tables ITT1 and ITTT]) ■ Confirming that our computational procedure can give reasonable 
results for bulk ZnO systems, we proceed to analyze ZnO single layers and nanotubes. 



III. ZINC OXIDE SINGLE LAYERS 

In this section, we discuss the elasticity and piezoelectricity of ZnO single layers. First, 
as an initial configuration, we take a single Zn-0 layer (shown in the left side of Fig. Q]) cut 
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from the wurtzite ZnO crystal, and then optimize its structure. As an stable configuration, 
we obtain a planar graphite-like structure as shown in the right side of Fig. ^ Here we 
use 6x6x1 Monkhorst-Pack k-points and two different groups of other input parameters 
for calculations. The final optimized results, listed in Table II VI are insensitive to the input 
parameters. The bond length of the planar structure is slightly smaller than that of ZnO 
crystals. Its binding energy is -8.246 eV/ZnO, higher than the value -8.947 eV/ZnO of ZnO 
crystals. 

Now we consider a wurtzite ZnO thin film with two polar surfaces (±0001 faces). Its 
energy is higher than that of bulk ZnO with the same number of atoms by E c = 4.0 J/m 2 
(the cleavage energy density^) because of the spontaneous polarization in c direction. Here 
we derive the the critical layer number based on our above result for planar ZnO single 
layers. Considering the influence of spontaneous polarization, we estimate the total binding 
energy E f = M(-8.947iV + E c n ) = M(-8.947iV + 2.21) eV for the wurtzite thin film with 
N Zn-0 single layers, where f2 is the area for a unit cell of the ZnO crystal in ab plane 
and M the number of unit cells in each layer. If we consider the thin film with N planar 
single layers and neglect the very weak interlayer attractions, the total binding energy is 
about -8.246MA eV. We find E f > -8.246MA for A < 4, which implies that a ZnO thin 
film less than N c = 4 Zn-0 layers prefers the planar graphite-like structure to the wurtzite 
structure. Here the estimated value 4 is the lower bound of critical number N c because we 
neglect the weak interlayer attractions. Note that the cleavage energy density may depend 
on the layer number and structural relaxation. Considering the structural relaxation, Meyer 
and Marj4i found that the the cleavage energy density varied from 3.0 J/m 2 to 3.4 J/m 2 , 
weakly dependent on the the layer number. Using these values, the lower bound of critical 
number N c is changed to 3. If we consider the interlayer attractions coming from the dipole- 
dipole interaction between planar ZnO layers because of their spontaneous polarizations 
(the ground state: spontaneous polarizations antiparallel in the neighbor layers), we should 
obtain iV c > 3, which is the exact meaning of "the lower bound of critical number being 3". 
Claeyssens et al. also obtained a critical thickness by DFT calculations.— They investigated 
the binding energy of the thin films with N = 4, 6, • • • , 24 single planar layers (layer distance 
2.4 A) and compared it with the binding energy of the thin films (wurtzite structure) with 
the same number of Zn-0 layers. They found that the ZnO thin film less than eighteen Zn-0 
layers prefers the planar graphite-like structure to the wurtzite structure. As is well known, 
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the widely used DFT packages cannot accurately describe the interlayer interactions . 43 ^ 44 
Thus the critical layer number N c = 18 obtained directly from DFT calculations might be 
inaccurate although it satisfies iV c > 3. With the development of DFT method, one may 
obtain the accurate critical layer number in the future when DFT packages can deal well 
with the interlayer interactions. 

Adopting the optimized structure, we can calculate the elastic and piezoelectric constants 
of ZnO single layers. From Fig. El we find that the ZnO single layer has the 3-fold rotation 
symmetry and a reflection symmetry respect to Xi-axis. Thus the stress-strain relation is 
expressed in the matrix form>^ 
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There are only two independent elastic constants: en and c\%. The piezoelectricity can be 
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where is the spontaneous polarization in the Xi-direction. There is only one independent 
piezoelectric coefficient en. The calculated values of the elastic and piezoelectric constants 
of ZnO single layers are listed in Table IS Using the elastic constants, we calculate the 
Young's modulus per atom by Y atom = (en — cf 2 /cn)/2 ~ 32.4 eV/atom, which is quite 
smaller than the value 57.7 eV/atom for a single graphite layer obtained in our previous 
work. 45 Additionally, we observe that the ZnO single layer (en = 0.48 e/A) has stronger 
piezoelectricity than the BN single layer (en = 0.23 e/A)4& 



IV. SINGLE- WALLED ZINC OXIDE NANOTUBES 

Generally speaking, if a material can exist in graphite-like structure, the corresponding 
nanotubes can be synthesized in the laboratory— A natural question is: how can we fabricate 
SWZONTs? We propose a possible way in this section. 

Without considering the two ends, an SWZONT can be thought of as a cylinder rolled 
up from a single sheet of ZnO layer such that two equivalent sites of the hexagonal lattice 
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coincide. To describe the SWZONT, we introduce several characteristic vectors in analogy 
with a single-walled carbon nanotube. As shown in Fig. El the chiral vector C^, which defines 
the relative location of two sites, is specified by a pair of integers (n, m) which is called the 
index of the SWZONT and relates to two unit vectors a x and a 2 of the hexagonal lattice 
(C/i = na.i + ma 2 ). The chiral angle 9 defines the angle between a x and CV For an (n,m) 
nanotube, 9 = arccos[(2n+m)/(2\/n 2 + m 2 + nm)\. The translational vector T corresponds 
to the first point in the hexagonal lattice through which the line normal to the chiral vector 
C/j passes. 

We calculate the binding energy of SWZONTs with different indexes which is shown in 
Fig. |21 From this figure we find that: (i) The binding energy for different SWZONTs is 
negative, which suggests these SWZONTs can exist in principle; (ii) The binding energy 
(Ef,) decreases with the increase of the radius (R) of SWZONTs and can be well fit by 

E b = -8.242 + 1.371/i? 2 (eV/ZnO). (5) 

Obviously, Ej, — > —8.242 eV/ZnO for R — ► oo. This value is quite close to the binding energy 
(-8.246 eV/ZnO) of a ZnO single planar layer. The term 1.371/i? 2 reflects the curvature 
effect of nanotubes. The classic shell theory also gives the same form, DQ/R 2 , for the 
curvature effect,^ where D is the rigidity of the ZnO single layer and = 8.9lA 2 is the area 
of the parallelogram generated by the unit vectors ai and a 2 . Thus we obtain D = 0.15 eV, 
which is quite smaller than the rigidity (1.17 eV)^ of a single graphite layer. The ZnO single 
layer is much softer than the graphite layer such that it should be more easily wrapped up 
into a nanotube. 

The solid-vapor phase process has been employed to successfully synthesize ZnO 
nanobelts, nanorings, and nanohelixes J^^ il(|) It might also be used to realize SWZONTs 
with carbon nanotubes as templates. As shown in Fig. the ZnO powder decomposes into 
gaseous Zn 2+ and O 2- at high temperature and low pressure. 49 After that, Ar flow carries 
Zn 2+ and O 2- to the low temperature zone at relative higher pressure where Zn 2+ and O 2- 
will deposit on the carbon nanotube array prepared carefully in advance. Controlling the 
low enough density of Zn 2+ and 2 ~ and the proper distance between carbon nanotubes, 
one may obtain thin enough ZnO films, for example less than 3 single layers, coating on 
the surfaces of carbon nanotubes. These films prefer the multi-walled tubes to the wurtzite 
structure because the layer number is smaller than 3. One may take out the SWZONTs 
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from the multi-walled tubes by analogy with the method in carbon nanotubes.— 

Now we discuss two points in the above process. The first point is the growth temper- 
ature. From the conditions for producing ZnO nanobelts , 7 i 8 i 9 i 10 the growth temperature 
of ZnO nanotubes is expected to be 300-500 °C. At this temperature, the highly purified 
carbon nanotubes are chemically inactive,— which ensures that they merely work as tem- 
plates but do not react with Zn 2+ or 2 ~. The other point is the critical radius of the 
synthesized SWZONTs. The carbon-ZnO multi-walled tubes can form in the above process 
if the attraction between nanotubes overcomes the curvature energy of ZnO nanotubes. For 
simplicity, we consider a carbon nanotube coated by an SWZONT. The attraction between 
them is approximately described by the Lennard- Jones potential. The layer distance is es- 
timated to be 3-4 A from the van der waals radii of carbon, oxygen, and zinc atoms. 52 The 
attraction energy (E at ) is about —19.5 meV/ZnO.— From E at + DQ/R 2 < 0, we obtain 
the smallest radius R c ~ 8.3 A synthesized in the above process. We should emphasize 
that the above process is merely possible method to produce the ZnO nanotubes in terms 
of the first-principles calculation and the previous experiments ^ 1 '? 11 ^ on the growth of ZnO 
nanobelts. Whether it works or not, experimental researches will gives an answer in the 
future. 

It is useful to discuss the piezoelectricity of SWZONTs. We need not to do additional 
DFT calculations considering the experience obtained by Sai et al. in the study of the 
piezoelectricity of single-walled BN nanotubes.™ They found that the piezoelectricity of BN 
nanotubes depends on the chirality. The physical origin is that an angle between local 
polarization of each unit cell and the tube axis depends on the chiral angle. Their key point 
is that the curvature of nanotubes has so small effect on the piezoelectricity that we can 
directly deduce their piezoelectric property from the piezoelectricity of single layers. Two 
deformation modes will result in the change of the polarization in the direction of tube axis: 
one is the tension or compression strain (ex) along the tube axis, another is the shear strain 
(7, the twist angle per length) induced by the torsion around the tube axis. The variation 
of the polarization can be expressed as 

SP T / L = e T e T + e t 7, (6) 

where L represents the perimeter of the nanotube. Following Refl46l one can derive the 
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coefficients 

er = en cos 39 (7) 

and 

e t = —en sin 39 (8) 

from the transformation between coordinate systems {xi,X2} and {C^, T} shown in Fig. [21 
From above three equations, we find that the piezoelectricity of SWZONTs depends on their 
chiral angles. For zigzag tubes (m = 0, 9 = 0°), only tension or compression strains induce 
polarizations, while only shear strains play roles for armchair tubes (n = m, 9 = 30°). 
Both strains work for achiral tubes. For a specific pure axial strain, the variation of the 
polarization decreases with the increase of chiral angle 9 from 0° to 30°. For pure shear 
strain induced by the torsion around the tube axis, the absolute variation of the polarization 
increases with the increase of chiral angle 9 from 0° to 30°. These effects are more evident 
in SWZONTs than in BN nanotubes because the ZnO single layer has a relative larger en 
than the BN single layer. 

V. CONCLUSION 

We have calculated the elastic constants and piezoelectric coefficients of ZnO crystals and 
single layers. We find that ZnO single layers are much more flexible than graphite single 
layers in the elasticity and stronger than boron nitride single layers in the piezoelectricity. 
SWZONTs are demonstrated to exist in principle and might be realized by the solid-vapor 
phase process with carbon nanotubes as templates. 

Finally, we would like to point out some potential applications of SWZONTs and related 
structures, (i) SWZONTs, same as other ZnO nanostructures, can be applied in optoelectric 
devices. The quality of SWZONT-based devices should be higher than others because of 
their very small thickness, (ii) The piezoelectricity of SWZONTs depends on their chirality. 
For most ZnO nanotubes except the zigzag type, twists around the tube axis will induce 
axial polarizations. This property makes SWZONTs possible to be used in the torsion mea- 
surement devices. For example, one may use SWZONTs to make very small Coulomb torsion 
balance, (iii) SWZONTs, except the armchair tubes, are polar tubes. 54 Water transport, 
ice formation, and biopolymer translocation in polar tubes should exhibit quite different 
characteristics from those in non-polar tubes. The latter has been fully studied, especially 
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for carbon nanotubes , 55 > 56 i 57 while the former has attracted no attention even for BN tubes 
up to now. (iv) The carbon-ZnO nanotubes formed in the solid-vapor phase process must 
have excellent mechanical and electrical properties: On the one hand, they have high enough 
rigidity and axial strength because of the contribution of carbon nanotubes. On the other 
hand, they exhibit strong piezoelectricity coming from the ZnO tubes. These structures 
may be employed to make electric nanotube motors^ where the axial voltage will be greatly 
reduced relative to carbon nanotube motors 58 because of the high piezoelectricity of ZnO 
tubes. 
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TABLE I: Structure parameters of ZnO crystals. 



Authors 


Methods 


a (A) 


c(A) 


c/a 


u 


Present 


DFT 


3.199 


5.167 


1.615 


0.379 


Karzel et al. a 


Experiment 


3.250 


5.204 


1.602 


0.382 


Wu et al. b 


DFT 


3.197 


5.166 


1.616 


0.380 


Catti et al. c 


Hartree-Fock 


3.286 


5.241 


1.595 


0.383 



"Reference SC . 
b Reference 23. 
c Reference |3J; ■ 



TABLE II: Elastic constants (relaxed-ions 


) of ZnO crystals 


in units of GPa. 




Authors 


Methods 


en 


C12 


Cl3 


C33 


C44 


Present 


DFPT 


218 


137 


121 


229 


38 


Wu et al. a 


DFPT 


226 


139 


123 


242 


40 


Catti et al. h 


Hartree-Fock 


246 


127 


105 


246 


56 


Kobiakov c 


Experiment 


207 


118 


106 


210 


45 


Azuhata et al. d 


Experiment 


190 


110 


90 


196 


39 



a Reference 23 ■ 
b Reference ■ 
c Reference S£ . 
^Reference 39. 
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TABLE III: Piezoelectric coefficients (relaxed-ions) of ZnO crystals in units of C/m 2 . 



Authors 



Methods 



e3i 



e33 



ei5 



Present 
Wu et al. a 
Catti et al. b 
Kobiakov c 



DFPT 
DFPT 

Hartree-Fock 
Experiment 



-0.65 
-0.67 
-0.55 
-0.62 



1.24 
1.28 
1.19 
0.96 



-0.54 
-0.53 
-0.46 
-0.37 



"Reference 
b Reference 
c Reference 



35. 

38. 



TABLE IV: Optimized structure of a ZnO single layer. Here L v , Bi, B a , and Eb represent the 
vacuum layer thickness, bond length, bond angle, and binding energy, respectively. 



L v (A) 



ecut (Ha) 



B i (A) 



B a 



E b (eV/ZnO) 



10.6 
21.2 



35 
60 



1.852 
1.852 



120° 
120° 



-8.246 
-8.247 



TABLE V: Elastic and piezoelectric constants (relaxed-ions) of ZnO single layers. Here L v repre- 
sents the vacuum layer thickness. 



L v (A) 


ecut (Ha) 


cu (eV/ZnO) 


cu (eV/ZnO) 


en (e/A) 


10.6 


35 


72.6 


24.0 


0.48 


21.2 


60 


72.7 


24.1 


0.48 



FIG. 1: (Color online) ZnO single layers. The left side is a single layer taken from the wurtzite 
ZnO crystal while the right one is the optimized structure which is a planar hexagonal lattice. The 
small and large balls represent zinc and oxygen atoms, respectively. 
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FIG. 2: (Color online) Unrolled honeycomb lattice of an SWZONT. By rolling up the sheet such 
that the atoms in the two ends of the vector coincide, a nanotube is formed. The vectors ai 
and a? are unit vectors of the hexagonal lattice. The translational vector T is perpendicular to 
and runs in the direction of the tube axis. 




R(A) 



FIG. 3: Binding energy (Ef,) and radius (R) of different SWZONTs. 



ZnO powder CNTarray 

t , — j decompose (9) deposit 

=C> 2 (g) = _ 



fl? 



High temperature zone 



Low temperature zone 



FIG. 4: (Color online) Schematic diagram of the vapor-solid phase process to synthesize SWZONTs 
with carbon nanotube (CNT) array as templates. 
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